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Abstract. We address the question of the “partition regularity” of the 
Pythagorean equation a^+b^ = ; in particular, can the natural numbers 

be assigned a 2-coloring, so that no Pythagorean triple (i.e., a solution 
to the equation) is monochromatic? We prove that the hypergraph of 
Pythagorean triples can contain no Steiner triple systems, a natural ob¬ 
struction to 2-colorability. Then, after transforming the question into one 
about 3-CNF satisfiability and applying some reductions, a SAT solver 
is used to find a 2-coloring for {1,..., 7664}. Work continues as we seek 
to improve the reductions and extend the computation. 
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1 Introduction 

One of the central problems in Ramsey Theory on the integers (see, for example, 
m) is the following question: given an equation - or set of equations - and a 
positive integer k, does there exist a coloring of N so that there are no monochro¬ 
matic solutions to the equation(s)? That is, suppose the family T of subsets of 
N is defined to be the set of solutions {xi ,..., a;„} to a system of equations in n 
variables. Does there exist a function c : N —>■ [A:] so that, for each F G T, there 
are Xi,Xj G F with c{xi) ^ c{xj)? A system of equations is called “partition 
regular” if, for any number of colors k, every fc-coloring of the integers admits a 
monochromatic solution. Perhaps the most famous examples of partition regu¬ 
larity are the equations x + y = z (Schur’s Theorem) and x + y = 2z (van der 
Waerden’s Theorem); Rado’s Theorem provides a vast generalization of these 
examples. 

A question that is surprisingly resistant to extant methods is the partition 
regularity of the “Pythagorean equation” x"^ + y"^ = z^; see, e.g., [5]. To date, it 
is not even known if it is possible to 2-color the naturals so that there are no 
monochromatic Pythagorean triples. However, there is a natural translation of 
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this question into that of solving a certain 3-CNF satisfiability problem: 

/\ {xi V Xj V Xk) A {-^Xi V -^Xj V -^Xk) ■ (1) 

i,j,k 

In this case, “true” and “false” assignments are the colors given to the indices 
of variables. 

In order to show that the Pythagorean triples are partition regular, it sufhces 
to demonstrate a subfamily of them which is not bipartite, i.e., any 2 -coloring 
induces a monochromatic triple. (The existence of such a family follows from the 
de Bruijn-Erdds Theorem, q.v. [3].) It is known that no two Pythagorean triples 
share two points, i.e., no leg and hypotenuse of one right triangle are the two 
legs of another integer right triangle. Indeed, the two equations x^ + y"^ = and 
yield two squares and z^) whose sum and difference are also 
squares. One can find an elementary proof that no such pair of squares exists in 
Sierpihski’s classic text [Ij. 

Therefore, in our search for nonbipartite subsystems of triples, we need only 
consider those in which each pair of triples intersect in at most one point. Such 
families are known as “linear 3-uniform hypergraphs” (in the graph theory com¬ 
munity) and “partial triple systems” or “packings” (in design theory). The small¬ 
est nonbipartite partial triple system is the ubiquitous Fano plane Fy, the projec¬ 
tive plane over F 2 with seven points and seven triples. Fy is also a “Steiner triple 
system,” meaning that each pair of points is contained in exactly one triple. In 
fact, it is no hard to see that every Steiner triple system is not bipartite, so it 
is reasonable to search for these classic obstructions to 2 -colorability among the 
Pythagorean triples. However, we show below that a broad class of triple systems 
including the family pyth of all Pythagorean triples contains no Steiner triple 
systems. Therefore, any search for nonbipartite partial triple systems in pyth 
necessarily must consider other systems than Steiner triple systems. 

In the next section, we introduce some notation and definitions. Section 
contains the proof that pyth, and all ordered triple systems with the “sum 
property”, contain no Steiner triple systems. In the following sections, we de¬ 
scribe the satisfiability solving methodology used to show that at least [7664] is 
2-colorable without any monochromatic Pythagorean triples. We conclude with 
some directions for further research. 

2 Preliminaries 

A triple system V. is a pair {V,E) = {¥{?{), E{'H)) consisting of a vertex set 
V and a family of unordered triples E C (^) . A partial triple system, aka a 3- 
uniform linear hypergraph, is a triple system in which each pair of distinct edges 
intersect in at most one vertex. A Steiner triple system is a partial triple system 
S so that, for each pair of vertices x and y, there is some z so that {x, y, z} is an 
edge of S. An ordered triple system % is a triple {V, E, <) = {V{H), E{H), <-u) 
consisting of a vertex set V, a family of unordered triples E C ( 3 ), and a total 
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ordering < on V. We say that an ordered triple system "H has the sum property 
if, whenever {a, 6 , c} and {a', b', c'} are two edges with respective maxima c and 

c', 

(a < a') A (6 < h') (c < c'). (2) 

Examples: 

1. Let E = Z, and, for each a ^ 6 , let e = {a,b,a + b) be an edge of E. We 
call this the Schur Triple System, and denote it by SCHUR. It is easy to see 
that SCHUR has the sum property. This shows immediately that the sum 
property does not imply finite colorability, since Schur’s Theorem says that 
any coloring y : N —>■ [c] with finitely many colors c admits a monochromatic 
e S SCHUR, i.e., |y(e)| = 1. Such a map is called a (weak) hypergraph color¬ 
ing. A map X : N —)■ [c] so that |x(e)| = 3 instead of just |x(e)| > 1 is known 
as a “strong coloring,” and corresponds exactly to a proper vertex coloring 
of the complement of the “leave,” i.e., the graph of all pairs contained in 
some triple e € E{'H). Some areas of mathematics call this the “shadow” 
graph or “1-skeleton” of TL, and denote it by the boundary operator dTL. 

2. It is clear that any order-preserving isomorphic image of a triple system 

with the sum property also has the sum property, and that any subhyper¬ 
graph of a triple system with the sum property does as well. For exam¬ 
ple, we define the Pythagorean Triple System PYTH by E(pyth) = N and 
if(PYTH) = {{a, b,c} : a'^ + b^ = c^}. Since one can embed pyth into SCHUR 
monotonically by the map , pyth has the sum property as well. It is 

a wide open problem to determine whether pyth has a weak coloring with 
finitely many colors. (It is even open whether it is strongly colorable.) As 
mentioned in the introduction, pyth is actually linear, i.e., no two edges 
intersect in more than one vertex. 

3. A special subsystem of pyth is the Primitive Pythagorean Triple System 
PRIM consisting of all Pythagorean triples which are relatively prime. That 
is, P = N and 

if(pRiM) = {{a, b, c} : a'^ + b^ = and gcd(a, b, c) = 1}. (3) 

It is easy to see that pyth is actually a union of dilates of prim by each d S N. 
However, prim is bipartite: the parity coloring n i— >■ n (mod 2) provides a 
2 -coloring. 

We define a special class of partial triple systems called “bicycles.” The k- 
bicycle has 2fc -I- 2 vertices and 2k edges. Its vertices are the elements of Z 2 fc and 
two “antipodes” a and b; its edges are all triples of the form {a, 2j, 2j -|- 1} and 
{b,2j —l,2j}, 0 < j < k. The 2-bicycle is also known as the Pasch configuration, 
or quadrilateral: the six-point partial triple system consisting of the edges abc, 
ade, bef, cdf. The 3-bicycle appears in the literature as the “hexagon” (e.g., 0): 
eight points {a, b, d, e, f, g, h, i\ with edges {afh, aei, adg, beh, bdi, bfg}. 

The following proposition follows immediately from well-known results in the 
theory of triple systems. (See, for example, 0.) For completeness, we give a short 
proof. 
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Fig. 1. A 3-bicycle and a 5-bicycle. 


Lemma 1. Ifv,w are vertices of nontrivial Steiner triple system S, then there 
is a k-bicycle for some k > 2 in S with antipodes v and w. 

Proof. Define the “link” Sx of a vertex x € 5 to be the set of pairs {a, b} so that 
{a, b, x} is an edge of S. For some z € S, {v, w, z} is a triple of S, so Sy consists 
of a perfect matching Mi on S — {v,w,z}, plus the edge {w,z}\ similarly, Sy, 
consists of a perfect matching M 2 on S — {v,w,z}, plus the edge {v,z}. The 
union of Mi and M 2 is composed of even-length cycles of length at least 4; any 
one of these forms a bicycle with antipodes v and w. 

We now define two weaker implicants of the sum property. We say that 
an ordered triple system has the upper sum property if, whenever {a,b,c} and 
{a,b',c'} are two edges with c and c' their respective maxima, 

(6 > 6') ^ (0 c'). (4) 

This clearly follows from the sum property by setting a = a'. We say that 
an ordered triple system has the lower sum property if, whenever {a, b, c} and 
{a',b',c} are two edges with c as both of their maxima, 

(a > a') ^ {b < b'). (5) 

To see that the lower sum property follows from the sum property, suppose that 
a > a' but b > b'. Then it follows that the maximal element of {a, b, c} is greater 
than the maximal element of {a',c}, whence c > c, a contradiction. We say 
that two pairs of integers {a, b} and {c, d} are “nesting” if a < c < d < 5. Then 
it is possible to restate the upper sum property as the fact that, for each vertex 
X, the subset of the fink graph Hx intersecting {y : y > x} is a non-nesting 
matching. The lower sum property may be similarly restated as the fact that. 
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for each vertex a;, the subset of the link graph T-Lx induced by {y : y < x} is a 
fully nested matching, i.e., for each two edges e and /, e is nested in / or vice 
versa. 

3 Sum Property implies No STS 

Proposition 1 If H has the upper sum property, and Q is a k-bicyele in H, 
then the maximal two points of Q are not its antipodes. 

Proof. Let 

Q = {Z 2 k U {a, b}, {aOl, a23, a45,...} U {612,634,656,...}) c 76 ( 6 ) 

be a A:-bicycle, and suppose a and 6 are the maximal two points of Q. We may 
assume without loss of generality that a > b. Then the maximal elements of all 
triples are a or 6 (depending on which antipode they contain). Therefore, since 
{a, 2 j, 2 j + 1 } n { 6 , 2j + 1 , 2j + 2 } = {2j + 1 }, 

(a > 6 ) ^ (2j > 2j + 2) (7) 

for each 0 < j < k. However, the quantities above are modulo 2k, whence the 
set of resulting inequalities is circular and therefore inconsistent. 

Proposition 2 If P has the lower sum property, and Q is a k-bicycle in %, 
then the maximal two points of Q are not its antipodes. 

Proof. Let 

Q = {1^2k U (a, 6 }, (oOl, a23, a45,...} U (612,634,656,...}) C 76 ( 8 ) 

be a fc-bicycle, and suppose a and 6 are the maximal two points of Q. We 
may assume without loss of generality that a > b. Then the maximal elements 
of all triples are a and 6 (depending on which antipode they contain). Since 
{a, 2 j, 2 j + 1} is an edge for each 0 < j < k, the k pairs {2j,2j + 1} are a 
matching, and they are linear ordered by nesting, by the lower sum property. 
Suppose, without loss of generality, that (0,1} is the outermost matching edge 
and 0 < 1, so that, for all x S Z 2 fc \ {0,1}, 0 < x < 1. The pairs { — 1,0} and 
{ 1 , 2 } are also nested, since they arise from the edges { 6 ,— 1 , 0 } and { 6 , 1 , 2 }. 
Hence, —1 < {1,2} < 0 or 1 < { — 1,0} < 2, and each of these possibilities 
contradicts 0 < 1 . 

Corollary 3 // 76 has the full, lower, or upper sum property, then it does not 
contain any Steiner triple system. In particular, PYTH contains no Steiner triple 
system. 

Proof. Suppose 76 contained some Steiner triple system T. Let a and 6 be the 
maximal elements of T. Then a and 6 are the antipodes of some bicycle, by 
Lemma[^ However, this contradicts Propositionsand/or 

Note that a triple system with the sum property ean contain a quadrilateral: 
for example, SCHUR contains {5,15, 20}, {5,8,13}, {7,8,15}, {7,13, 20}. 
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4 The Computational Approach 

We now describe the methodology used to certify the 2-colorability of the induced 
hypergraph of pyth on the vertex set [7664], 


4.1 Listing the Triples up to N 

To make a list of all Pythagorean triples having only elements less than or equal 
to some upper bound N, we use a function based on Dickson’s Method, from 
Wikipedia [7]. Our function for generating triples uses the equations = 2st 
for r, s, and t G N, and x = r + s,y = r + t and z = r + s + t, with x'^ +y'^ = . 

With these equations, we step through all combinations of s and t so that r G N 
and z < N. The Python code is short. 

import numpy as np 
def dicksonGenTriples(uBound): 
data = [] 

sSteps = range(1,uBound) 
for s in sSteps: 

tSteps = range(s,uBound-s) 
for t in tSteps: 

r = np.sqrt(2*s*t) 
if r == np.round(r): 
r = int(r) 
z = r + s + t 
if z >= uBound: 

break 
X = r + s 
y = r + t 

if [x,y,z] not in data: 
data.append{[x,y,z]) 

with open('triples-dickson.json', 'w') as f: 

json.dump(data, f) 
return data 


4.2 Writing the CNF 

We remapped the integers involved in the Pythagorean triples to consecutive 
integers starting with 1 since, some SAT solvers require it. The bijection for 
remapping is represented in the Python code by a dictionary. 

{" 3 " : 1 , " 4 ": 2 , " 5 ": 3 , " 6 ": 4 , " 8 ": 5 , " 10 ": 6 } 

The following is a CNF file for the triples with integers up to 10. These are [3,4,5] 
and [6,8,10]. 
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c 10 

p cnf 6 4 
12 3 0 
-1 -2 -3 0 
4 5 6 0 
-4 -5 -6 0 

A full specification for CNF (Conjunctive Normal Form) is given at [8]. 

4.3 Choosing a Solver 

We experimented with glueSplit_clasp, Dimetheus, SWDiASBY, PeneLoPe, Lin- 
geling, and Treengeling, all from http://www.satcompetition.org/ in 
2014. Treengeling is a parallel solver that runs on some or all of the cores of 
a single machine, and seemed to be the fastest for a CNF corresponding to the 
set of triples up to about 7000. These runs completed within a few hours on each 
of the many solvers. We did not emphasize finding the best solver. 

4.4 Reducing the Size of the Input to the Solver 

A natural goal of SAT solving is to ease the burden on the solver by finding 
parts of the hypergraph, corresponding to the solver’s original CNF, which can 
be safely removed without changing the satisfiability of the new, reduced CNF. 
We have had to rely on these reductions for upper bounds above 7620. 

Our only success to date at reducing the size of the input to the solver in 
this way, is “pendant removal”. By pendant, we mean any edge (triple) of the 
3-uniform hypergraph which has at least one vertex not contained in any other 
edges. Once a triple system with a pendant removed is properly 2-colored, the 
pendant’s vertex, or vertices, of degree one can then be assigned any color that 
is not assigned to some other vertex of its edge (of which there is always at least 
one). Pendants are always removable and can be removed iteratively until no 
more pendants remain. Once these pendants have been removed from a graph 
and the remaining graph has been fed to the solver and colored, it is easy to add 
back the removed pendant edges and 2-color them as they are added back. The 
edges are added back in the reverse order that they were removed. 

5 Parallelizing Across Multiple Computer Cluster Nodes 

5.1 The Overall Strategy 

We use a strategy to parallelize across multiple nodes of a cluster similar to one 
used in a paper by Biere, Heule, Kullmann, and Wieringa [9]. In particular, we 
convert the problem into 2”^ many potentially easier problems, by making 2™ 
copies of the original CNF file and appending a distinct list of extra clauses to 
each one. We run these new CNF files on up to 2™ cluster nodes; one node per 
CNF. As soon as any node completes, if that CNF is satisfiable, we have our 
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solution and can force all the other nodes to terminate before completing. The 
various nodes do not communicate amongst themselves and run independently 
from one another. 

To create the extra clauses for each of the 2™ copies of the original CNF 
file, we select m special vertices (integers), each coming from some Pythagorean 
triple which was used to generate the original CNF. For each special vertex, i, 
we make a new clause 

{^Xi V -^Xi V -^Xt) (9) 

or 

{xi V Xj V Xi). (10) 

This gives us 2™ different lists of m extra clauses to add to each of the 2™ 
copies of the original CNF file. These extra m clauses provide the solver with 
predetermined truth assignments for specific vertices. Either one of these 2™ runs 
of the SAT solver, on the slightly varying CNF files, will complete and provide 
a valid coloring, or all the runs will complete and each certify unsatisfiability - 
given sufficient computational time. 

5.2 An Example CNF File Modified for Parallelization 

We present one of the 4 CNF files for the list of triples [[3,4,5],[6,8,10]], split on 
2 vertices; 3 and 4, which are remapped to 1 and 2. The line which begins with 
p is also updated by adding m = 2 to the last integer in the original CNF file 
on that line. The original CNF had a 4 in that spot. Note especially the last 2 
clauses. 

c 10 

p cnf 6 6 
12 3 0 
-1 -2 -3 0 
4 5 6 0 
-4 -5 -6 0 
1110 
-2 -2 -2 0 

5.3 Choosing Vertices to Split On 

It would be desirable to choose a list of m special vertices so that each vertex 
in the list is relatively independent from the others. By independence, we mean 
that any assignment of truth values to the variables corresponding to special 
vertices leads to a similar number of (partial) satisfying assignments. We have 
attempted to gauge independence by measuring the time needed for each of the 
2™ truth value combinations to run using Treengeling. The minimum distance 
in the hypergraph of all Pythagorean triples, up to an upper bound, N, between 
any two of the special vertices is a reasonable proxy for their independence. 
Under current technical limitations, we have found m = 2 to run the fastest. 
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The BFS Method of Finding Distant Vertices A promising method of 
choosing some vi and V 2 which are distant from one another in the hypergraph 
of Pythagorean triples is to use what we call the breadth-first search (BFS) 
method. By way of illustration, we present Figure and Table We start from 
a seed triple to, and use the method to find a tn that is furthest from tg- Once a 



Fig. 2. A graph of some hypothetical list of triples where any 2 triples share an edge 
if they intersect. For example to = [3,4,5] and ti = [5,12,13] share an edge in a graph 
of this type, as the 2 triples share a vertex. 


Table 1. The triples of each BFS level corresponding to the graph in Figurej^ 

level 0123456 7 
to tl t2 to t4 ts to til 
tr tg tg tio 


triple appears in a level, it will not appear again in the list. When there is more 
than one path from the seed to some other triple, the listing in Table shows 
only the shortest path. Thus any two triples which are in two distant levels from 
one another, necessarily have a high minimum graph distance between them. 

In the Pythagorean Triple System, to find two vertices that are relatively far 
apart, we choose a seed triple with minimum constituent integers as small as 
possible. In the example of Figure]^ we run the BFS method on a list of triples 
that has had its pendants removed. The plot gives an indication of the proximity 
of two vertices to each other and for the connectedness, loosely defined, of each 
of the two vertices to the rest of the hypergraph. 
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Level 


Fig. 3. The bar plot shows the size of each BFS level for the list of pendant-removed 
triples with all vertices less than or equal to 7664. 


Gauging the Independence of a List of Special Vertices As an illustration 
of a methodology to gauge the independence of a set of vertices, we choose a 
list of 4 randomly selected special vertices and uniformly at random remove 
10 % of the triples in the list of triples with vertices less than or equal to some 
upper bound around say, 7621. This speeds up the SAT solver computations 
from about a day to about 100 seconds. Each run uses 1 node with 12 cores. 
Next we record the run time for 10 runs each, of the 16 truth value combinations, 
for the 4 special vertices. Each of the 10 runs of the SAT solver uses a different 
reduced list of triples. There are a total of 10 reduced triple lists for the whole 
process. See Figure]^ For each of the 16 truth value assignments, we compute 
an empirical mean run time; then, the sample variance of these 16 numbers is 
then computed. For a set of chosen special vertices, a low variance indicates 
high independence, since the solver worked approximately equally hard across 
the different possible truth values. 


5.4 Running a Node Pool on Clusters 

We have devised a method for simultaneously running a fixed number of nodes, 
each running an instance of the solver on 12 cores. We are allowed to use up to 
about 10 of these nodes at a time, after waiting in the queue for some amount 
of time. See a detailed description of our clusters at the University of South 
Carolina, Maxwell and Planck, at 
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100 



True / False Combination Number 


Fig. 4. Rnn times of 10 runs of each of the 16 possible truth values of the 4 special 
vertices. The special vertices were chosen uniformly at random. 


http://WWW.sc.edu/about/offices_and_divisions/division_of_ 
information_technology/rci/hpc_resources/index.php 

We use the two clusters interchangeably. Their specifications are as follows. 

— Maxwell 

• Hardware 

* 40 GL390 Nodes: 12 cores per node, Intel Xeon 2.4 GHz, 24 GB 
RAM 

* 6 SL250: 16 cores per node, Intel Xeon 2.60GHz, 32 GB RAM 

* 1 DL380 Headnode: 12 core, 48GB RAM 

* 24TB attached storage 

• Interconnect 

* QDR Infiniband 

• Software 

* Centos 

* HP CMU cluster management utility 

* OpenMPI 

* Torque/Maui scheduler 

— Planck 

• Hardware 

* 20 SL250 Nodes: 12 cores per node, Intel Xeon 2.8 GHz, 24 GB RAM 

* 15 X 3 NVIDIA M1060 with 240 cores each, 3x3 

* NVIDIA M2070 with 448 cores each 

* 11 TB attached storage 
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• Interconnect 

* QDR Infiniband 

• Software 

* Centos 

* HP CMU cluster management utility 

* OpenMPI 

* Torque/Maui scheduler 

Our parallelization starts from an outer shell script which runs a sequential 
Python program that creates a fixed number of inner shell scripts, each running 
1 instance of Treengeling on one full node. The Python program, which is started 
by the outer shell script, uses system calls to detect the completions of each of 
the inner shell scripts by detecting the existence of the output file for each inner 
shell script. When an inner shell script completes, a new inner shell script is fed 
to the queue until all the desired instances of Treengeling have been submitted, 
or until one instance of Treengeling returns a valid coloring. 


5.5 Choosing the Number of Special Vertices 

We present some evidence that splitting approximately 4 ways using 2 special 
vertices is fastest, given current technical constraints. See Figure which plots 
run time against the number of nodes at which the problem was split. 



Fig. 5. Some results of scalability testing. 
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6 Results 

6.1 Performance of Treengeling 

Treengeling can find solutions for all the Pythagorean triples with integers up 
to 7620 in about 21 hours. If the list of triples has been pendant removed, 
Treengeling can find a solution for all triples up to 7650 in about 33 hours. 


6.2 Freedom of Coloring and Searching for Patterns 

We have looked for some patterns in the large valid colorings found, but have 
not found much. There seems to be a lot of freedom in how a coloring can be 
achieved, which seems to be a reason why we have not found any consistent 
patterns in the colorings. 


6.3 A Visualization of the Coloring up to 7664 

The coloring in Figure was made by removing the pendants from the list of 
triples with integers up to 7664, running the solver on the result, and adding 
back the triples, one at a time, in the reverse order that they were removed and 
2-coloring them as they were added back. The final coloring was checked to see 
that every triple with all integers less than or equal to 7664 is dichromatic. 

Integers which do not occur in any triple, with all constituent integers of the 
triple being less than or equal to 7664, are white. The solver colors every integer 
which it receives from any triple as either false or true. While some integers 
might be colorable as either false or true in a particular coloring, the solver will 
assign some color to each of these integers along with every other integer that it 
receives. 

7 Future Work 

7.1 Further Reducing the Input 

We are running a sequential Python function on the clusters that systematically 
checks all possible linear 3-uniform hypergraphs on n vertices to see if any fur¬ 
ther reductions to pyth can be made. We have not yet found any “removable” 
subgraphs that are not pendants, from the triple list up to 7700. We would like 
to parallelize this sequential operation using Python and MPI, or C and MPI if 
need be. 


Acknowledgements Thanks to Ron Graham, Bill Kay, and Christopher Poirel 
for helpful discussions and ideas in the development of the present work. Paul 
Sagona has been helpful with the high performance computing resources at the 
University of South Carolina. 
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Fig. 6. A coloring of the hypergraph of pyth induced by the vertex set [7664]. False is 
grey. True is black. Integers not appearing in the input to the solver are white. The 
integer positions climb vertically by columns from the lower left to the upper right. 
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